Valence, spin, and orbital state of the Co ions in the one-dimensional Ca3Co206: an 
x-ray absorption and magnetic circular dichroism study 
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We have investigated the valence, spin, and orbital state of the Co ions in the one-dimensional 
cobaltate Ca3Co206 using x-ray absorption and x-ray magnetic circular dichroism at the Co-L2,3 
edges. The Co ions at both the octahedral Cooct and trigonal Cotrig sites are found to be in a 
3+ state. From the analysis of the dichroism we established a low-spin state for the Cooct and a 
high-spin state with an anomalously large orbital moment of l.lfiB at the Cof^g ions. This large 
orbital moment along the c-axis chain and the unusually large magnetocrystalline anisotropy can 
be traced back to the double occupancy of the d2 orbital in trigonal crystal field. 

PACS numbers: 71.27.-l-a, 71.70.Ej, 75.20.Hr, 75.30.Gw 
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The one-dimensional Ca3Co206 has attracted great in- 
terest in recent years due to the observation of the stair- 
step jumps in the magnetization at regular intervals of 
the applied magnetic field i'^ii^i^i^i^ii^'^ The rhombohc- 
dral structure of this compound consists of [Co206]oo 
chains running along the c axis of the hexagonal unit 
cell.^° In each chain, CoOg octahedra alternate with 
CoOg trigonal prisms. The magnetism is Ising-like and 
directed along the Co chains with large magnetic mo- 
ments of 4.8 /iB per formula unit.^^:i3. The intra-chain 
coupling is ferromagnetic with a spin-freeze at Tgp = 7 K 
and the chains couple antiferromagneticly with a Neel 
temperature of Tjv = 25 Kfii 

Based on density-functional-theory calculations, Vidya 
et al. claimed a low-spin (LS) Co^+ and a high -spm 
(HS) Coj+g state for Ca3Co206.^^ However, other ex- 
perimental and theoretical works have proposed a Co^+ 
valence state at both the octahedral and trigonal Co 
sites, with the Cooct LS {S = 0) and the Cotrig HS 
(5 = 2) state J^ii^ii^ii^ The Ising character of the mag- 
netism is also subject of discussion. Dai and Whangbo 
et al. found from their band structure calculations that 
the spin-orbit-inactive do orbital lies lowest of all Cotrig 
crystal-field levels, -i^iiP and had to invoke excited states 
in their attempt to explain the Ising magnetism. Wu et 
al-r^ on the other hand, calculated that it is the spin- 
orbit-active d2 orbital which lies lowest, giving a very 
different picture for the Ising magnetism. 

In this paper we have applied soft-x-ray absorption 
spectroscopy (XAS) and magnetic circular dichroism 
(XMCD) at the Co-L2,3 edges to resolve the Co valence, 
spin and orbital state issue in Ca3Co206. We have also 
carried out detailed configuration-interaction cluster cal- 



culations to analyze the spectra. We found that the Co 
ions are all in the 3-1- state, that the Cooct are non- 
magnetic, and that the Co^+g carry a 1.7^b orbital mo- 
ment. We clarify the orbital occupation issue and the 
origin of the Ising magnetism. 

The Ca3Co206 single crystals were grown by heating a 
mixture of Ca3Co409 and K2CO3 in a weight ratio of 1:7 
at 950 °C for 50 h in an alumina crucible in air. The cool- 
ing was performed in two steps, first down to 930 °C at a 
rate of 10 °C/h and then down to room temperature at 
100 °C/h4The Co-L2,3 XAS spectra of Ca3Co206 and of 
CoO and EUC0O3 (Ref. [10) as references were recorded 
at the Dragon beamline of the National Synchrotron Ra- 
diation Research Center (NSRRC) in Taiwan with an 
energy resolution of 0.3 eV. The first sharp peak at 777 
eV of photon energy at the C0-L3 edge of CoO was used 
for energy calibration, which enabled us to achieve better 
than 0.05 eV accuracy in the relative energy alignment. 
The XMCD spectra were collected at the ID08 beamline 
of the European Synchrotron Radiation Facility (ESRF) 
in Grenoble with a resolution of 0.25 eV and a degree 
of circular polarization close to 100% in a magnetic field 
of 5.5 Tesla with the sample kept at approximately 15- 
20 K, using a dedicated superconducting magnet setup 
from Oxford Instruments. The Poynting vector of the 
photons and the magnetic field were both parallel to the 
c axis. The single-domain Ca3Co206 crystal used for the 
XMCD experiment was needle-shaped with a dimension 
of 0.2 X 0.2 X 10 mm"^ for a x b x c. Clean sample areas 
were obtained by cleaving the crystals in situ in chambers 
with base pressures in the low 10"^" mbar range. The Co 
^2,3 spectra were recorded using the total-electron-yield 
method (TEY). 
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FIG. 1: Co-L2,3 XAS spectra of Ca3Co206, CoO, and 
EuCoOs. 



Fig. [T] shows the Co-L2,3 XAS spectrum of Ca3Co206 
together with that of CoO and EuCoOa. CoO serves here 
as a 2+ reference and EuCoOa as a LS 3+ referencei^ 
One can see first of all that the 2+ spectrum (CoO) con- 
tains peaks which are 2 or more eV lower in energy than 
the main peak of the 3+ spectrum (EuCoOa). It is well 
known that x-ray absorption spectra at the transition- 
metal ^2,3 edges are highly sensitive to the valence state. 
An increase of the valence of the metal ion by one re- 
sults in a shift of the L2,3 XAS spectra to higher en- 
ergies by 1 eV or more.^^ This shift is due to a final 
state effect in the x-ray absorption process. The energy 
difference between a 3(i" {3(f for Co^+) and a 3cr'~^ 
(3d^ for Co^+) configuration is AE = E{2p^3(P-'^ 
2p^W)-E{2p^W 2p53d"+i) « Upd-Udd « 1-2 eV, 
where Udd is the Coulomb repulsion energy between two 
3d electrons and Upd the one between a 3d electron and 
the 2p core hole. 

One can observe from Fig. [T] that the Ca3Co206 spec- 
trum has no features at the low-energy side, which oth- 
erwise could have indicated the presence of Co^+ species 
like in CoO. Instead, the spectrum has a much closer re- 
semblance to that of Co'^+, like in EUC0O3. This means 
that one can safely rule out the Co^+/Co^+ scenario.— 
In other words, the XAS experiment reveals unambigu- 
ously that both the Cooct and the Cotrig ions are in the 
3-1- valence state. This result supports the analysis of 
the Co 2p core-level x-ray photoemission spectra^S. and 
band-structure calculation s^ "^'--'. 

In order to resolve the spin-state issue, we now resort 
to XMCD. The top part (a) of Fig. [2] depicts 00-^2,3 
XAS spectra taken with circularly polarized light with 
the photon spin parallel (red solid line, and antipar- 
allel (black dashed line, aligned to the magnetic field. 
The quantization axis z has been chosen to be parallel 
to the c axis, which is the easy magnetization axis"^ . One 
can clearly observe large differences between the spectra 
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FIG. 2: (Color online) (a) Measured soft-x-ray absorption 
spectra with parallel (^i"*", red solid line) and antiparallel {fJ,~ , 
black dashed line) alignment between photon spin and mag- 
netic field, together with the difi'erence spectrum (/i^ — fi~ , 
blue dash-dotted); (b) Simulated sum spectra assuming a dou- 
bly occupied d2 orbital for the Cotrig and low-spin (LS) Cooct 
ions; (c) and (d) Contribution of the Cotrig and Cooct ions to 
the simulated sum spectra. 



using these two alignments. The difference spectrum, i.e. 
the XMCD spectrum, is also shown (blue dash-dotted 
line, fj,'^ — fj,~). Using the well-known XMCD sum rule 
developed by Thole et al.r^ 



where iVe is the number of the Co d electrons, we can 
extract directly the orbital (Lz) contribution to the mag- 
netic moment without the need to do detailed modeling. 
Assuming an average Co 3d occupation number of about 
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Nf. = 6.5 electrons (estimated for a HS Co^"*" oxide^^ and 
also to be justified below by cluster calculations) we find 
a value of 1.2/iB for the L^, which is a very large num- 
ber indicating that the ground state is strongly spin-orbit 
active. 

To extract more detailed information concerning the 
spin and orbital states from the Co L2,3 XAS spectra, we 
have carried out simulations for the XMCD spectra us- 
ing the well-proved configuration-interaction (CI) cluster 
model.— i^^i^ The method uses for each Co site a CoOg 
cluster which includes the full atomic multiplet theory 
and the local effects of the solid. It accounts for the 
intra-atomic 3c?-3d and 2p-3d Coulomb interactions, the 
atomic 2p and Zd spin-orbit couplings, the O 2p-Co id 
hybridization, and the proper local crystal-field parame- 
ters. In the configuration-interaction cluster calculation 
we describe the ground state by the configurations 

= aeld*^) + a^\d'L) + as\d^l/) + agld^L^), (2) 

where L denotes a ligand hole and J2^=6 = 1 '^'''^^ The 
transition-metal electron occupation is then given by 

Ne^6al + 7a^ + 8al + 9al. (3) 

The simulations have been carried out using the program 
XTLS 8.3.24 

In octahedral symmetry the 3d orbitals split up in the 
well-known Cg and t2g levels; the splitting is given by 
lODq. In a trigonal prismatic environment, however, it is 
found that the — is degenerate with the xy^ and the 
yz with zx orbitah ^^'^^ In the presence of the spin-orbit 
coupling, it is then better to use the complex orbitals c?Oj 
dijd-i, and di/d^i. Band structure calculations indi- 
cate that the di/d^i band is split off from the do a-nd 
^2/^-2 bands by about 1 eV, while the do f^nd d2/d-2 
bands are almost degeneratei ^^i^^ Critical for the mag- 
netism and for the line shape of the simulated spectra are 
the crystal field parameter lODq for the CooctOg cluster 
and the crystal field parameter for the splitting between 
the nearly degenerate do and d2 orbitals of the CotrigOg 
cluster. lODq needs to be critically tuned since this is 
set to determine whether Cooct-Og cluster is in the LS 
or HS state^^. The trigonal prism crystal field has also 
to be fine tuned since it determines whether the d2 or the 
do lies lowest, and thus the magnitude of the orbital mo- 
ment and strength of the magnetocrystalline anisotropy 
as we will show below. Tuning of these parameters will 
be done by establishing which of the simulated spectra 
reproduce the experimentally observed ones. 

As a starting point, we make the assumption that the 
Cooct ion is in the LS state, based on the observation 
that the average Cooct-0 bond length of 1.916 A in this 
materiali^ is shorter than the 1.925 Ain LaCoOa at 5 K, 
which is known to be LS-^^. With such a short length, 
the Co ion is subjected to a large enough lODq, sufficient 
to stabilize the non-mag netic LS state^^i^o for Coltt 
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in Ca3Co206. With this starting point, the magnetism 
and the XMCD signal have to originate from the Cotrig 
ions. This is quite plausible since with the extremely 
large Cotrig-0 average bond length of 2.062 A,-'^^ which 
is much larger than 1.961 A for LaCoOs at 1000 K^O 
can expect that the crystal field is small enough to sta- 
bilize the HS state. 2- Based on the observation that the 
orbital contribution to the magnetic moment is extremely 
large, we take the d2 ansatz for the Cotrig and not the 
do. In the CI calculation, the parameters for the multi- 
pole part of the Coulomb interactions were given by the 
Hartree-Fock values, while the monopole parts (Udd, Upd) 
were estimated from photoemission experiments on Co'^"'" 
materials2i. The one-electron parameters such as the O 
2p-Co 3d charge-transfer energies and integrals were es- 
timated from band-structure results. iiii^ii^i The charge- 
transfer energy is given by A = E{d^) — E{d'^L) = 1.5 eV, 
the d-d Coulomb repulsion by Udd = 5.5 eV and of p- 
d of the excited Co by Upd = 7.0 eV; the Slater inte- 
grals have been reduced to 80% of their Hartree-Fock 
value. An exchange field of 7?ex = 3 meV has been 
used. For the Cotrig ions the ionic crystal splittings 
are AEl°i'j = 0.9 eV and AEl°Z = 0.05 eV taken 

"l/"2 "0/"2 

from band-structure calculation's,; the hybridization is 
V^]'y^ = 1.88 eV, V^]"^^ = 1.28 eV, and V^^^ = 1.25 eV. 
For Cooct ions pda — —1.44 eV and pdn — 0.63 eV was 
used. With this set of parameters we have found a LS- 
HS transition for Cooct at IQDq — 0.65 eV. Here we 
used lODq = 0.8 eV, based on band-structure calcula- 
tion, which is the same value as for EuCoOs known as a 
LS Co3+ oxide.20 

The results of these simulations are given by curves (b) 
in Fig. [21 together with a break-down into the separate 
contributions of the Cotrig (curves c) and Cooct (curves 
d) ions. One can clearly observe that the simulations 
(curves b) reproduce the experimental spectra quite well. 
All major and minor peaks in the individual experimen- 
tal XAS (/i"*", fi~) and XMCD (/i"*" — ^~) spectra (curves 
a) are present in the simulations (curves b). It is almost 
needless to remark that the entire simulated XMCD sig- 
nal is coming from the HS Cotrig ions (curves c), since 
we started with a nonmagnetic LS Cooct (curves d). In 
the simulation, we find Sz = 1.8/xb and = l-T/is, giv- 
ing a total magnetic moment of 2Sz + Lz — 5.3fiB per 
formula unit. Due to strong hybridization between Co 
3d and O 2p, the 3d occupation numbers of the Cotrig 
ions and the Cooct are 6.3 and 6.8, respectively, giv- 
ing on average 6.5 electrons as used for the sum rules. 
The calculated total moment from the simulation com- 
pares reasonably well with the 4.8/is value from direct 
magnetization measurements.-'-^ Yet, the simulated Lz 
value {1.7 fis) is appreciably larger than the one derived 
from the experiment using the XMCD sum rule (1.2/is). 
However, looking more closely at the simulated and ex- 
perimental XMCD spectra, we can clearly see that the 
XMCD spectra have very similar line shapes and that the 
distinction is mainly in the amplitude, i.e. a matter of 
scaling. This indicates that the discrepancy in the Lz de- 
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FIG. 3: (Color online) (a) Measured soft-x-ray magnetic cir- 
cular dichroism spectrum (XMCD, /i"*" — (b) Simulated 
XMCD spectrum assuming a doubly occupied do orbital for 
the Cotiig and low-spin (LS) Cooct ions. 



termination might be caused by the fact that the sample 
is not fully magnetized in our experiment. According to 
magnetization measurements one can only achieve 90% 
of the saturation magnetization at 5.5 T^i^ In addition, 
slight misalignment of the sample together with the high 
magnetocrystalline anisotropy may account for some fur- 
ther reduction of the experimental value. 

Having established that the (i2/LS scenario for the 
Cotrig/Cooct ions explains well the experimental spectra, 
we now investigate the sensitivity of our analysis. For 
this we change the ansatz for the Cotrig ion: we now put 
the do orbital to be energetically lower than the d2'?- The 
result is shown in Fig. O in which the simulated XMCD 
spectrum (b) is compared with the experimental one (a) . 
One can unambiguously recognize large discrepancies in 
the line shapes, not only in the L3 region (777-785 eV), 
where the simulated XMCD signal has much less ampli- 
tude, but also in the L2 (792-797 eV), where now an 
XMCD signal is calculated while it is practically absent 
in the experiment. Using the XMCD sum rulc;^ we can 
relate these discrepancies also directly to the fact that the 
do ansatz essentially does not carry an orbital moment. 

The success of the cluster method for the analysis of 
both the high-energy spectroscopies and the ground-state 
magnetic moments provides confidence for its use to in- 
vestigate the magnetocrystalline anisotropy in this mate- 
rial. We have calculated the total energy of the CotrigOg 
cluster as a function of the exchange field H^x, directed 
either along the z axis (the c axis or the one-dimensional 
chain) or along the x axis (perpendicular to the c axis). 
The results for the ^2 ansatz are plotted in the left panel 
of Fig. m It is evident that the lowest state gains en- 
ergy when Hex is along z and very little to nothing for 
flex along X. This demonstrates directly that the mag- 
netocrystalline anisotropy should be large and directed 
along the z, consistent with the experimental observa- 
tion that the magnetization parallel to c is almost sixfold 
of the one perpendicular c (at 12 K and 5.5 T)2, and with 
the Ising nature of the magnetism. We note that the 
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FIG. 4: (Color online) Total-energy level diagram for the 
Cot rig ion as a function of the exchange field H^x along the 
z direction (c axis, solid black lines) and along the x direc- 
tion (perpendicular to the c axis, dashed blue lines), with (left 
panel) the d2 and (right panel) the do orbital doubly occupied. 



rate twice that of Hex (along z), meaning that Sz is close 
to which in turn is consistent with the HS {S = 2) 
nature of the Cotrig ion. 

We have also analyzed the magnetocrystalline 
anisotropy of the CotrigOg cluster using the do ansatz. 
The right panel of Fig. |4] reveals that the total energy 
decreases for both directions of H^x, but now with the 
important distinction that the energy for H^x along x 
is always lower than for Hex along z. This implies that 
the easy-magnetization axis should be perpendicular to 
the z axis, which is not consistent with the experimental 
facts. Also this contradiction thus effectively falsifies the 
do ansatzJ^^ 

To conclude, using soft x-ray absorption spectroscopy 
and the magnetic circular dichroism therein at the Co- 
^2,3 edges we have experimentally determined that both 
Co ions in Ca3Co206 are 3-I-, with the Cotrig ion in the 
high-spin state and the Cooct in the nonmagnetic state. 
The Cotrig ion carries an anomalously large orbital mo- 
ment of 1.7 fj,B which we have been able to relate to the 
double occupation of the d2 orbital. In addition, the de- 
tailed analysis of the spectral line shapes together with 
that of the magnetocrystalline anisotropy firmly estab- 
lishes that the d2 orbital lies lowest in energy^^ and not 
the do.— 1^ This in turn also demonstrates that a proper 
incorporation of the spin-orbit interaction is required for 
the ab-initio calculation of the delicate electronic struc- 
ture of this material. 
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